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Exploring the magnetized cosmic web through low frequency radio 
emission 
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ABSTRACT 

Recent improvements in the capabilities of low frequency radio telescopes provide a unique 
opportunity to study thermal and non-thermal properties of the cosmic web. We argue that the 
diffuse, polarized emission from giant radio relics traces structure formation shock waves and 
illuminates the large-scale magnetic field. To show this, we model the population of shock- 
accelerated relativistic electrons in high-resolution cosmological simulations of galaxy clus- 
ters and calculate the resulting radio synchrotron emission. We find that individual shock 
waves correspond to localized peaks in the radio surface brightness map which enables us to 
measure Mach numbers for these shocks. We show that the luminosities and number counts 
of the relics strongly depend on the magnetic field properties, the cluster mass and dynami- 
cal state. By suitably combining different cluster data, including Faraday rotation measures, 
we are able to constrain some macroscopic parameters of the plasma at the structure forma- 
tion shocks, such as models of turbulence. We also predict upper limits for the properties of 
the warm-hot intergalactic medium, such as its temperature and density. We predict that the 
current generation of radio telescopes (LOFAR, GMRT, MWA, LWA) have the potential to 
discover a substantially larger sample of radio relics, with multiple relics expected for each 
violently merging cluster. Future experiments (SKA) should enable us to further probe the 
macroscopic parameters of plasma physics in clusters. 

Key words: magnetic fields, cosmic rays, radiation mechanisms: non-thermal, elementary 
particles, galaxies: clusters: general 



understood. The magnetic fields influence the physics of the plasma 
in several important ways. They couple the collisionless charged 



. 1 INTRODUCTION AND KEY QUESTIONS through magneto-hydrodynamic turbulent interactions with cosmic 

■ rays leading to second order Fermi acceleration. They illuminate 

r > The plasma within and between galaxies is magnetized. Despite ... , i-.-i.i_t i 

Vh 1 distant cosmic ray electron populations by enabling synchrotron 
rri many observational efforts to measure galactic and intergalactic . . , . ,, ... .. , " , 

. . emission and tell us indirectly about violent high-energy astrophys- 

magnetic fields, their properties and origins are not currently well . , . „ . . " , _, 

• ical processes such as formation shock waves or y-ray bursts. The 

magnetic fields in spiral galaxies are highly regular, showing align- 
ment with the spiral arms. They are believed to arise from weak 
particles to a single but complex fluid through the Lorentz force, , " ,■ , ,_r • , 

seed fields amplified by dynamo processes, driven by differential 
and trace dynamical processes in the Universe. Magnetic pres- .... r™ , n ,, ,,, , , , 

J r or rotation in galactic disks. The seed fields could have been produced 

sure and tension mediate forces and provide the plasma with ad- , " r . „ . , , . . „ 

„ _ ..,„,, by many sources, ranging from stellar winds and jets of active 

ditional macroscopic degrees of freedom in terms of Alfvemc . . , . . ...... , , _ ■ , , 

galactic nuclei, to plasma instabilities and battery effects in shock 
and magnetosonic waves. They cause the turbufent cascade to be- " . . . .. . . . .... 

c J waves, in ionization fronts, and in neutral gas-plasma interactions, 

come anisotropic towards smaller scales and suppress transport 

processes such as heat conduction and cosmic ray diffusion across 



More hypothetical ideas for the seed field origins invoke primordial 

generation in early Universe processes, such as phase transitions 
the mean magnetic field. They are essential for accelerating cos- , . , „ . „ .. T , , . , , 

° J ° during the epoch of inflation. In order to understand more about 

mic rays by providing macroscopic scattering agents which en- * , . . , , , , 

„ . magneto-genesis, we need to study the least processed plasma pos- 

ables diffusive shock acceleration (first order Fermi process) and ... ... .... , ,. . __,.... 

sible that still shows some primordial memory. This points us to 

the magnetized plasma in intergalactic space, in particular to the 

plasma in galaxy clusters. There, magnetic fields show a smaller 
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ing by turbulent gas flows, driven by galaxy cluster mergers, and 
the orbits of the member galaxies. For an overview on the present 
observational and theoretical kno wledg e the reader is pointed 
by Reesl 



to the review articles 



1987) 



1 19931); iKronberd dl994h: iBecketalJ 



Wielebinski & Krause 
19961); iKulsrud J 1999' 



Beckl(l200lh:lGrasso & Rubinsteid ( l200il) : ICarilli & Tavlod < l20oL 
Widrow (2002). This work aims at closing a gap between theoreti 



cally motivated phenomenological models of large scale magnetic 
fields and actual observational non-thermal phenomena associated 
with them. 

Diffuse radio synchrotron e mission has alread y been observed 
in more than 50 galaxy clusters dFerrari et alj|2008l) . The emission 
is associated with the entire intra-cluster medium (ICM). The syn- 
chrotron emission process demonstrates the presence of highly rel- 
ativistic electrons (cosmic ray electrons, CRe) with a Lorentz factor 
typically up to 7 ~ 10 4 and magnetic fields within the ICM. The 
diffuse radio emission can be classified into two categories: radio 
halos and radio relics. Giant radio halos are centrally located, trace 
the thermal emission and show no sign of polarization, while radio 
relics are located at the periphery of clusters, are polarized and are 
elongated in appearance. There exist a number of classes of radio 
objects that have bee n referred to over the years as "radio relics" 
(Kem pner et al.l2004l and references therein). Two of these are as- 
sociated with extinct or dying active galactic nuclei (AGN). These 
either host a synchrotron cooling radio plasma from a past AGN 
outburst that created the radio lobes or are revived "radio ghosts" 
where an aged r adio relic has been re-energized by a merger or an 
accretion shock jEnfflin & Gopal-Krishnj|200lh . 

The focus of this paper is on a third type of radio relic 
emiss ion, sometimes referred to as radio "gischt'Q ( Kempner et al. 
2004), that shows diffuse emission on scales up to 1 Mpc. Diffu- 
sive shock acceleration at structure formation shocks can energize 
a primary p opulation of relati v istic electrons that emit synchrotron 
radiation (EnBlinetal. 1 1 19981 : 1 Minia ti et al. I l200lh in a magnetic 
field that can be amplified by the post-shock turbulence. Promi- 
nent e xamples for this class o f radio relics c an be seen in Abell 
3667 dRottgeringetalJ Il997l) Abell 2256 dBridle & Fomalontl 
19761 : iMasson & Maverll 19781 : iBridle et al.|[l979l: iRottgering et al.1 
1994 : IClarke & Enfflinl 120061) . Abell 3376~ilBaechi et alj |2006|) . 

recently, Abell 2255 jPizzo et alj|20081) and Abell 521 



and, more 
dGiacintucci 



.centiy. 

etalJ l": 



2008). All galaxy clusters with observed radio 



relic emission are merging or show signs of ongoing dynamical ac- 
tivity, but not all dynamically active galaxy clusters are observed to 
have relics. This raises the question of whether diffuse radio emis- 
sion is a property of a special subset of clusters, or a universal prop- 
erty, with many relics too faint to be seen by current telescopes. 

From CMB measurements we kn ow that the Universe is com- 
posed of 4.6% baryonic matter, e.g. Ko matsu et al. (2008). How- 
ever, when observing the local U niverse (z < 1), we can account for 
fewer than half of these baryons (Fukugita 2004; Danf orth & Shulll 
2005). This is known as the missing baryon problem. The cur- 
rent cosmological paradigm of large scale structure formation pro- 
vides a solution to the missing baryon problem. As the Universe 
evolves, large scale structure grows from small density perturba- 
tions imprinted during an earlier epoch. In the hierarchical scenario 
of structure formation, structure grows from small to large scales, 
with baryons flowing in a filamentary web with clusters at the in- 



1 The name "gischt" derives from a German word for the crest on top of 
waves that are breaking at the shore thus resembling the radio emission of 



freshly injected electrons by formation shocks. 



terstices. The temperature of baryons deviates from adiabatic cool- 
ing associated with the Hubble expansion by increasing multiple 
times in discrete steps - always corresponding to a passage through 
a structure formation shock. Before they are shock-heated to the 
virial temperatures kT ~ 1 - 10 keV of galaxy groups and clusters, 
where they can be observed through their thermal bremsstrahlung 
emission, they are predicted to reside in the warm-hot intergalac- 
tic medium (WHIM). Temperatures in the WHIM are in the range 
of 10 5 K < T < 1 7 K Sellsten et al ll 19981: ICen & Ostrikeill 19991 : 
iDave et alj|200ll : iFurlanetto & Loebll2004l : iKang et al.ll2005l) . We 
will investigate whether it is possible for diffuse radio emission as- 
sociated with these formation shocks to be used as a tracer of the 
WHIM boundaries and so indirectly observe the WHIM. 

A larger sample size of these diffuse radio sources is re- 
quired to delve deep into the details of the non-thermal pro- 
cesses working within galaxy clusters. However, the combina- 
tion of low surface brightness, diffuseness and small dynamical 
range in the sensitivity of current radio telescopes makes the de- 
tection of this particular radio emission difficult. With the cur- 
rent capabilities of th e Giant Meter Radio Telescope (GMRT, 
Ananthakri shnan 1 19951) and the imminent arrivals of the Low Fre- 
quency ARray (L OFAR. iRQttgeringl 12 003), the Murchison Wide- 
field Array ( MWA.lMorales et alj|200^7 and the Long Wavelength 

Array (LWA. lKassim et alj|2005l) . and eventually the construction 

I — — II 1* 

of the Square Kilometre Array (SKA, Keshet et al. 2004b), power- 
ful low frequency radio telescopes are positioned to further increase 
our understanding of diffuse radio emission and give us insight into 
the following important topics: 

• the strength and coherence scale of magnetic fields on scales 
of galaxy clusters, 

• the process of diffusive shock acceleration of electrons, 

• the existence and properties of the WHIM, 

• the exploration of observables beyond the thermal cluster 
emission which are sensitive to the dynamical state of the cluster. 

In the course of this work we will consider how radio relic 
emission can shed light on each of these topics. To do this, we 
adopt a simplified model for the shock-accelerated population of 
electrons. The key figure illustrating these considerations is shown 
in Fig. [T] In our simulations, we can visualize properties of struc- 
ture formation shocks that are triggered by a recent merger of a 
large galaxy cluster (M ^ 10 15 /i~' M s ) and dissipate the associ- 
ated gravitational energy. In the left panel, the shock Mach num- 
bers, weighted by the energy dissipation rate, are encoded by the 
colour, while the brightness displays the logarithm of the dissipa- 
tion rate. This shows that most of the energy is dissipated in weak 
flow shocks internal to the cluster, while the shock waves become 
strongest and steepen as they break at the shallower peripheral po- 
tentials of the clusters and within filaments. The right panel shows a 
three-color image of of energy dissipation rate at shocks (shown in 
red and yellow) and radio synchrotron emission at 150 MHz from 
shock-accelerate d relativistic electrons (blue and contours, mod- 
eled according to lPfrommer et alj|2008l) . This "radio gischt" emis- 
sion traces structure formation shock waves and highlights the in- 
termittent nature of mass accretion, in particular along the filament 
extending from the cluster center to the upper left of the image 
and for the giant radio relic to the lower-right of the cluster. This 
radio emission illuminates magnetic fields that are amplified by 
magneto-hydrodynamic instabilities that are associated with these 
shock waves. The paper has been structured accordingly: in Sect. [2] 
we outline our methodology; we describe our results and discuss 
them in Sects. [3]and|4j and present our conclusions in Sect. [5] 
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Figure 1. Structure formation shocks triggered by a recent merger of a large galaxy cluster (M - 10 i5 /j~' M Q ) dissipate the associated gravitational energy. 
Left: the Mach number of shocks weighted by the energy dissipation rate is shown by the colour (while the brightness displays the logarithm of the dissipation 
rate). Right: three-color image of energy dissipation rate at shocks (shown with a color scale ranging from black over red to yellow) and radio synchrotron 
emission at 150 MHz from shock-accelerated relativistic electrons (blue and contours with levels starting at 7 X 10~ 4 mJy arcmin~ 2 and increasing with a 
factor of 15, respectively). This "radio gischt" emission traces structure formation shock waves, highlights the intermittent nature of mass accretion in galaxy 
clusters, and illuminates magnetic fields that are amplified by turbulence that can be excited by these shock waves. 



2 METHODOLOGY 

We briefly summarize our procedure. We model the synchrotron 
emission by calculating the primary shock-accelerated electron 
population using a scheme that is based on the thermal leakage 
model - a model that has been developed in the context of diffu - 
sive shock acceleration at supernova remnants fellison etaljl98ll) . 
We use a simple parametrization for the magnetic field. This lets 
us quickly scan the observationally allowed parameter space asso- 
ciated with the mostly-unknown spatial distribution of shocks on 
cluster scales and beyond. In the post-processing, we search for 
spatially correlated synchrotron emission from formation shocks, 
which represent our simulated radio relics and study the properties 
of these relics in the clusters in our sample. Our aim is to under- 
stand how radio observables can be used to reconstruct the physi- 
cal properties of radio relics, which trace structure formation and 
large scale magnetic fields. In the subsequent sections there will be 
a detailed description of our simulations and modelling. 

2.1 Adopted cosmology and simulated cluster sample 

Our work is based on high resolution smoothed particle hydrody- 
namics (SPH) simulations of galaxy clusters (mi nimum gas mass 
resolu t ion ~ 8 x IO'/i 'Mq for more details, cf. IPfrommer et aU 
I2007L I 2008) using the 'zoomed initial conditions' technique 
(Kat z&WhrteT l993) that were selected from a low resolution dark 
matter only simulation jYoshida et al.ll200~ll ) with a box size of 479 
h~ l Mpc. They were carried out using a mo dified version of the 
massively parallel tree SPH code GADGET-2 dSpringell200l) . The 
simulations of the galaxy clusters were performed in a "Concor- 
dance" cosmology model, ACDM with cosmological parameters 
of: n m = fi D M + fib = 0.3, Cl b = 0.039, fi A = 0.7, h = 0.7, n s = 1 



and erg = 0.9. Here, Q m refers to the total matter density in units 
of the critical density today, p crit = 3H^/(&nG). £\ and Q. A de- 
note the densities of baryons and the cosmological constant at the 
present day. The Hubble constant at the present day is parametrized 
as H = 100 h km s _1 Mpc _1 , while n s denotes the spectral index of 
the primordial power-spectrum, and erg is the rms linear mass fluc- 
tuation within a sphere of radius 8 /?~'Mpc extrapolated to z = 0. 

The simulations include a prescription for radiative cooling, 
star formation, supernova feedback and a formalism for detect- 
ing structure formation shocks and measuring the as sociated shock 
strengths, i.e. the Mach numbers jPfrommer et alj I2006I) . Radia- 
tive cooling was computed assuming an optically thin gas of pri- 
mordial composition (mass-fraction of X a = 0.76 for hydrogen 
and 1 - Xu = .24 for he l ium) i n collisional ionisation equilib- 
rium, following iKatz et alj d 1996T) . We also included heating by 
a photo-ionising, time-dependent, uniform ultra violet (UV) back- 
groun d expected from a population of quasars faaardt &~M adau 
1996), which reionises the Universe at z - 6. Star formation 
is treated using the h ybrid multiphase model for t he interstellar 
medium introduced by Springel & Hernquist d2003ri . In short, the 
ISM is pictured as a two-phase fluid consisting of cold clouds that 
are embedded at pressure equilibrium in an ambient hot medium. 

The cluster sample is displayed in Table [T] From this sam- 
ple the cluster g72a was chosen for detailed analysis of the prop- 
erties of radio relics since it is a relatively large (with a mass 
M a 10 15 M o ) post-merging cluster, similar to the Coma cluster. 
Additionally, it hosts the brightest radio relic in the entire sample. 
This relic resembles already observed ones. 
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Figure 2. Surface brightness emission map for radio relics in the simulated cluster g72a Left to right: the emission from two sets of relic finder parameters, and 
the total primary radio emission at 150 MHz. Our relic finder groups SPH particles using a friends-of-friends algorithm; we additionally require these particles 
to exceed a emissivity threshold (ET). The differences in the images illustrate the dependence on ET (left panel: observable parameters with ET = 10 -43 h? 
erg s~' Hz~' ster~' cm -3 , central panel: theoretical parameters ET = 10~ 55 h 3 erg s _I Hz 1 ster~' cm -3 ). The central map only lacks surface brightness at 
the level of 10~ 3 compared to the total primary emission. Both relic emission maps to the left contain more than 99% of the total flux from the total primary 
emission map. 



Table 1. Cluster sample considered in this paper. 



cluster 
name 


dynamical 
state" 




nb 

n 2aa 
[hr x Mpc] 


KJ 200 

[keV] 


g8a 


CC 


1.8 x 10 15 


2.0 


13.1 


gla 


cc 


1.3 x 10 15 


1.8 


10.6 


g72a 


PostM 


1.1 x 10 15 


1.7 


9.4 


g51 


CC 


1.1 x 10 15 


1.7 


9.4 


gib 


M 


3.7 x 10 14 


1.2 


4.7 


g72b 


M 


1.5 x 1() 14 


0.87 


2.4 


glc 


M 


1.4 x 10 14 


0.84 


2.3 


gld 


M 


9.2 x 10 13 


0.73 


1.7 


g676 


CC 


8.8 x 10 13 


0.72 


1.7 


g914 


cc 


8.5 x 10 13 


0.71 


1.6 



" The dynamical state has been classified through a combined criterion in- 
voking a merger tree study and the visual inspection of the X-ray bright- 
ness maps. The labels for the clusters are M-merger, PostM-post merger 
(slightly elongated X-ray contours, weak cool core region developing), CC- 
cool core cluster with extended cooling region (smooth X-ray profile). b The 
virial mass and radius are related by Ma(z) = |?r Ap cr j t (z)R 3 , where A = 
200 denotes a multiple of the critical overdensity p c rit(z) = 3//(z) 2 /(8;rG). 
c The virial temperature is defined by kT& = GM^/im p /(2R^), where /i 
denotes the mean molecular weight. 



2.2 Realization of magnetic fields 

Current SPH implementations that are capable of following the 
magneto-hydrodynamics (MHD) of the gas are presently still 
fraught with numerical and physical difficulties, in particular 
when following dissipative gas physics dDolagetal.lll999l 120051 : 
IPrice & Monaghar] |2004l, |2005). Hence we apply a parametriza- 
tion in the post-processing of our completed simulations in order 
to determine the strength and morphology of the magnetic field 
jPfrommer82 008). Secondly, the parametrization approach provides 
us with the advantage of exploring the parameter space of our mag- 
netic field description more efficiently, since we are not required to 
re-simulate when we alter the ab-initio unknown magnetic field pa- 



rameters. We have chosen a simple scaling model for the magnetic 
field of 

(1) 



Our independent model parameters are the magnetic decline a B , 
and the magnetic core energy density e Bo . The thermal energy den- 
sity e th is measured in units of its central energy density e tho = 
3P(h /2, which we calculate by fitting a modified /J-model (Eqn.[2]l 
to the radial pressure profiles of our clusters. We first remove the 
over-cooled core (see Sect. I2.4.21 , 



P(X) 



th„ 



1+1^ 

r r 



(2) 



We found that our modified /?-model provides a better fit to the 
pressure profiles than the us ually adopted spherically symmet ric 
King profiles, i.e. a /J-model dCavaliere & Fusco-Fernianall978h . 

This parametrization (Eqn. [Q was mo tivated by non-radiative 
SPH MHD simulations jDolag et alj [l999) and radiative adaptive 
mesh refinement MHD simulations jDubois & Teyssier 200a) of 
the formation of galaxy clusters in a cosmological setting. Rather 
than applying a scaling with the gas density as those simulations 
suggest, we chose the energy density of the thermal gas. Current 
cosmological radiative simulations (that do not include feedback 
from AGN) over-cool the centres of clusters, giving an overpro- 
duction of stars, enhanced central gas densities, and lower central 
temperatures than are seen in X-ray observations. In contrast, the 
thermal energy density of the gas is well-behaved in simulations. 
Observationally, the parametrization (Eqn.QJ is con sistent with sta- 
tistica l studies of Faraday rotation measure maps JVogt & Enfilin 
2005). Theoretically, the growth of magnetic field strength is de- 
termined through turbulent dynamo processes that will saturate at 
a field strengt h determined by the strength of the magnetic back- 
reaction (e.g. Subramanian 2003; Schekochihi n & Cowlevlr200r3) 



2 Recently, Sijacki et al. 12008) found that including cosmic rays from 
AGN in SPH simulations can solve the over-cooling problem while pro- 
viding excellent agreement of the gas fraction and the inner temperature 
profile. 



© 2008 RAS, MNRAS 000,[T]|T8] 



Exploring the magnetized cosmic web 5 



Table 2. Magnetic field parameters in various combinations. We define our 
standard magnetic field parameters to be ag = 0.5, Bo = 5\iG and v = 
150MHz. These parameters are used throughout the paper unless otherwise 
stated. 





BoInG] 


magnetic decline 


core magnetic field strength 


0.3 


2.5 


0.5 


5.0 


0.7 


10.0 


0.9 





and is typically a fraction of the turbulent energy density. The tur- 
bulent energy density should be related to the thermal energy den- 
sity, thus motivating our model theoretically. The parameter e Bo is 
constrained by past measurements of magnetic fields within clus- 
ters and is c hosen such that Bp = [e Bn 8^] l/2 to be on the order 
of a f ew txG jGovoni et al1l200fj ; lTavlor et alJl2007l ; lGuidetti et all 
2008). The parameters explored in our model are shown in Table|2] 

To predict the polarization angle and Faraday rotation mea- 
sure in our simulations, we nee d to model th e magnetic mor- 
phology of the ICM. We follow Tribble (1991) in order to cre- 
ate individual components of the magnetic vector field that obey 
a given power spectrum. The details of the magnetic field struc- 
ture within the ICM are still unknown. There have been mea- 
surements of magnetic correlations from Faraday rotation measure 
(RM) maps which are however limited by the finite window size 
of radio lobes and hence only constrain the spectrum on smaller 
scales. These measurements suggest that the fields are tangled 
with a Kolmogorov/Oboukhov-type power spectru m for coherence 
lengths of approximately 10 kpc scales and smaller (Vogt&EnBlin 
2003, 20051; iGuidetti et alj2 008). It has been argued that shallower 
magnetic field pow er spectra allow for longer coherence le ngths on 
the order 100 kpc jMurgia et al.ll2004l ; lGovoni et alj|2006h . On the 
other hand, a Fourier analysis of XMM-Newton X-ray data reveals 
the presence of a scale-invariant pressure fluctuation spectrum in 
the range between 40 and 90 kpc and is found to be well described 
by a projected Kolmo gorov/Oboukhov-type turbulence spectrum 
( Schu ecker et al]|2004l) . Assuming that the growth of the magnetic 
field strength is determined through turbulent dynamo processes 
suggests a similar spectrum of magnetic and hydrodynamic turbu- 
lence (e.g. Subramanian 20030 ■ 

We model the components of the magnetic field, S ; , as ran- 
dom Gaussian fields. We use a Kolmogorov power spectrum on 
scales smaller than the coherence length, and a flat (white-noise) 
power spectrum on larger scales. All three components of the mag- 
netic field are treated independently, which ensures that the final 
distribution of B(r) has random phases. After mapping our SPH 
Lagrangian energy density distribution of the thermal gas onto a 
3D grid (cf. Appendix Eqn. lAll ). these realizations of the magnetic 
field are then scaled such that the magnetic energy density obeys 
our assumed scaling given by Eqn.[T] To ensure V ■ B = 0, we ap- 
ply a divergen ce cleaning procedure to our fields in Fourier space 
(Bals aral 19981) : 



3 / fc'lc \ 

Bm = 2 k ■- y- \Bj(k). 



Applying this procedure to our Gaussian random field removes a 
third of the magnetic energy. Thus, we re-normalize B to conserve 
the magnetic energy. 



2.3 Cosmic ray electrons and synchrotron emission 

Collisionless cluster shocks are able to accelerate ions and electrons 
in the high-energy tail of their Maxwellian distributi on functions 
through diffusive shock acceleration (for reviews see jDrurvll 19831 : 
iBlandford & Eichlerlll987l ; iMalkov & O'C Drurvll200ll) . Neglect- 
ing non-linear shock acceleration and cosmic ray modified shock 
structure, the process of diffusive shock acceleration uniquely de- 
termines the spectrum of the freshly injected relativistic electron 
population in the post-shock region that cools and finally dimin- 
ishes as a result of loss processes. The radio synchrotron emitting 
electron population cools on such a short time scale T sync < 10 8 yrs 
(compared to the very long dynamical time scale T dyn ~ 1 Gyr) that 
we can describe this by instantaneous cooling. In this approxima- 
tion, there is no steady-state electron population and we would have 
to convert the energy from the electrons to inverse Compton (IC) 
and synchrotron radiation. Instead, we introduce a virtual electron 
population that lives in the SPH-broadened shock volume only; this 
is defined to be the volume where energy dissipation takes place. 
Within this volume, which is co-moving with the shock, we can use 
the steady-state solution for the distribution function of relativistic 
electrons and we assume no relativistic electrons in the post-shock 
volume, where no energy dissipation occurs. Thus, the cooled CR 
electron equilibrium spectrum can be derived from balancing the 
shock injection with the IC/synchrotron cooling: above a GeV it is 
given by 



fe(E) = C e E~ 



C, oc 



P 



SB + £ph 



(4) 



Here, a e = ar in j + 1 is the spectral index of the equilibrium electron 
spectrum and e pll denotes the photon energy density, taken to be 
that of CMB pho tons. A mo re detailed description of our approach 
can be found in Pfrommer (2008|). The synchrotron emissivity j v 
for a power-law spectrum of CRes scales as 



h • 



C C S" 



(5) 



where a v = (a e - l)/2. A line-of-sight summation of j v yields the 
radio surface brightness, S v . The surface brightness are provided in 
units of h = 0.7 to simplify comparison with observations. 



2.4 Finding radio relics 

In our search for radio relics in the simulated clusters, we have 
modified a friends-of-friends (FOF) dGeller & Huchrall 19831) algo- 
rithm so that it groups together connected radio synchrotron emis- 
sion in 3D. This relic finder works in a manner similar to a FOF 
finder except we have introduced the additional criterion of an 
emission threshold which the SPH gas particle are required to ex- 
ceed before being assigned into a group. Thus, our algorithm de- 
pends on three internal parameters which determine the groups of 
particles that are designated relics: the linking length, emissivity 
threshold, and minimum number of particles (Fig. |2j. The linking 
length (LL) is the parameter which controls the maximum distance 
(d max ) between two particles that can still be considered neighbours, 



(3) 



^DM Pcrit . 



(6) 
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Figure 3. Left two panels: simulated GMRT map of the cluster g72a at z=0.05 (similar to the Coma cluster), with the reconstructed beam in the bottom-left 
corner. The left panel is the 'dirty' map and the right panel is the 'clean' map, where we removed the brightest relic found in the dirty map, mimicking the 
cleaning procedure of radio maps. Right panel: surface emission map of the relics of g72a using our calibrated set of observable parameters shown with 
a logarithmic colour scale. Note that we reproduce the relic emission that can principally be detected by GMRT after applying a cleaning procedure to the 
compact and diffuse radio emitting sources. 



where pdm = PctU^dm is the mass density of dark matter and 
(M DM jj is the average mass of our dark matter particles. The link- 
ing length and the emission threshold parameters have degenerate 
effects on the resulting groups of particles. Through inspection, 
we have chosen to fix the linking length value of 0.2 resulting in 
<4iax = 50 kpc and vary the emission threshold. A minimum par- 
ticle value of 32 regulates possible SPH shot noise and allows for 
smaller structures to be included in the relic catalogue. The final 
parameter in our relic finder, which we have chosen to vary is the 
emission threshold. We compute the synchrotron emissivity of all 
the particles using Eqn.[5]and compare it to the emission threshold. 
The sets of grouped particles for each of our clusters are our relic 
catalogues that form the basis of our study. 

2.4.1 Determination of the emission threshold 

We tailored the calibration of the emission threshold in the relic 
finder to two cases: firstly, so that we would find relics observable 
by GMRT and LOFAR; secondly, so that we could study the com- 
plete picture that might be achievable with future radio telescopes 
such as SKA, by pushing the emission threshold back to the limit 
of our simulations. For this procedure, we left the magnetic field 
parameters unchanged. (We used our standard magnetic field pa- 
rameters cf. Table[2]) 

We simulate the visibilities and maps from the relics that 
GMRT would observe. We use the GMRT primary beam and an- 
tenna positions projected against the zenith {i.e. ignoring the z- 
component of the antenna positions). For simplicity, we approxi- 
mate the continuous UV tracks for each baseline by circles in the 
UV plane with measurements 18° apart. We observe cluster g72a 
surface emission at a redshift of 0.05 (cf. Fig.|3j and make a 'dirty' 
radio map, the Fourier transform of the visibilities. We used an inte- 
gration time of 2.5 minutes with a sensitivity of 0.2 mjy Vhr/beam 
at the frequency v = 150 MHz in simulating the visibilities. 

To approximate GMRT's dynamical range, we modelled a 
simple cleaning procedure by removing the brightest relic in the 

3 Note that the quantity in the brackets is equivalent to the ratio of (M(,)/pb 
except that the baryonic phase consists of gas and stars. 



'dirty' map from the total surface emission map and re-simulated 
the GMRT detections - resulting in our 'cleaned' radio map. We 
compared these maps to surface brightness maps of different relic 
catalogues where we varied our emission threshold (cf. Fig.[Ul. The 
emission threshold which reproduced the simulated images most 
accurately was 10~ 43 h 3 erg s -1 Hz -1 ster -1 cm~ 3 with the linking 
length and minimum number of particles already fixed. These pa- 
rameters are referred throughout this paper as observable parame- 
ters. 

The choice of a second emission threshold is related to the 
peak of the emissivity distribution function which is determined 
by the mass resolution of the SPH particles in the simulations. We 
found the peak emissivity to be at j v = 10 -55 h 3 erg s~' Hz - ' ster~' 
citT 3 . We find that changing j v by six orders of magnitude does 
not change the number of relics in a significant way (cf. Appendix 
IBV making the difference between peak emissivity and the observ- 
able parameters reasonable. Together with the linking length and 
minimum number of particle parameters stated above, these param- 
eters are referred to as theoretical parameters. Since the emissivity 
scales with frequency, the emission threshold must scale with fre- 
quency as well. The emission threshold (ET) scaling is fixed at our 
reference frequency vo = 150 MHz, 

ET = ET ^ , (7) 

and ETq adopts the values quoted above for both observational and 
theoretical parameters. 

In summary, the observable parameters were chosen to pro- 
duce relic catalogues resembling the ones obtainable from current 
or near future observations, whereas the theoretical parameters 
lead to hypothetical catalogues only obtainable with a perfect 3- 
d tomography of the medium which may find application with the 
future radio interferometer SKA (cf. Table[3}- 

2.4.2 Removal of galaxy contamination and cool core 

Our radiative s i mula tions model star formation 
jHernq uist & Springel 120031) which leads to the formation of 
galaxies. When applying the criteria described in Sect. 12.41 these 
galaxies appear as false radio relic candidates since they are dense, 
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Figure 4. Luminosity functions with observable parameters (left) and theoretical parameters (right) of our relic finder. The top and bottom panels show how 
the magnetic decline ag and the magnetic core energy density Bo impact the luminosity functions, respectively. The choice of magnetic field parameters has 
a large impact on the shapes of the luminosity function which is a consequence of the inhomogeneous nature of virializing processes in cosmic structure 
formation. Smaller o-g values, corresponding to a shallower magnetic decline, produce many more brighter relics compared to larger values of ag, which 
produce slightly more less-luminous relics. Increasing Bo results in a greater number of more-luminous relics. Notice the loss of the less-luminous relics 
from the observable to theoretical parameters. This is a result of the brightest relic swallowing up smaller relics due to a decrease of the emission threshold 
parameter. 



Table 3. Parameters chosen for the relic finder. 



parameter 


linking 


minimum number 


emission threshold 


name 


length 


of particles 


[erg s~' Hz" 1 ster -1 cirT 3 ] 


observable 


0.2 


32 


io- 55 


theoretical 


0.2 


32 


io-« 



Our relic finder also picks up the over-cooled centres of galaxy 
clusters in our simulations, contaminating the radio emission. Since 
the candidate relic in the over-cooled center may be physically con- 
nected to other true relics, it cannot be removed by simply discard- 
ing the closest relic candidate to the center. We apply a very conser- 
vative cut in radius of r = 40 kpc and neglect the weak dependence 
on cluster mass and dynamical state. We note that smaller clusters 
(M < 5 x 10 14 M o ), in particular those with dynamical activity, tend 
to have slightly smaller cooling regions. 



compact and have enough emissivity per particle to be selected 
by the relic finder. To select against those objects, we impose 
further constraints on the SPH particles that are grouped together 
and require them to have a zero fraction of neutral hydrogen 
and to be below a very conservative threshold of number density 

"thres = «SFt/32 = 0.004 CITT 3 . 
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3 RESULTS 

3.1 Probing the intra cluster magnetic fields 

In this section, we investigate how sensitive different radio syn- 
chrotron observables are with respect to the properties of the large 
scale magnetic field. 



3.1.1 Luminosity functions 

The magnetic fields within our simulations are parametrized by a 
simple scaling relation. For each cluster, we compute radio lumi- 
nosity functions to aid in differentiating between different magnetic 
field parametrizations by employing the dependency of synchrotron 
emissivity on the magnetic field of the ICM. Our luminosity func- 
tions are distribution functions of the total luminosity per relic (J v ), 
where 



Jv ^ jv.a ^ Jy,a- 



(8) 



The units of J va are erg s -1 Hz 1 ster -1 , M a and p„ are the SPH gas 
particle mass and density respectively for the set of SPH particles 
within the relic labelled by a. 

The number of relics seen depends on the magnetic field 
parametrization. In Fig. [4] we show how the luminosity functions 
depend on B and a B . As expected, we find more and brighter radio 
relics for higher values of So. However, rather than simply scal- 
ing the luminosity function to higher relic emissivities for larger B 
(assuming a fixed slope a B ), we find that their shapes change. This 
is a consequence of the inhomogeneity of the virializing cosmic 
structure formation waves that are illuminated by the synchrotron 
emitting electrons. The effect of varying B is analogous to the wa- 
ter level within a very inhomogeneous landscape that corresponds 
to the strength of the virializing shock waves. This level can adopt 
different values depending on the magnetic realization such that 
the resulting synchrotron emitting objects end up single connected 
or disjoint. So, one could consider using Minkowski functionals to 
characterize the different relics. 

The trend for a B is the opposite: higher values of a B lead to a 
lower number of radio relics. The parameter a B represents the slope 
of the magnetic scaling (£ tho > £ th ), and a steeper slope will result in 
the magnetic field strength falling off faster with radius. We expect 
the effect of both So and a B on the luminosity function to be a 
generic effect for all the simulated clusters, since j v cc B a " +l /(e B + 
£cmb) (cf. Eqns.[4]and[5} and in the peripheral cluster regions where 
e b < £ cmb> we obtain j v oc S" ,+1 . The luminosity functions alone 
are not sufficient to fully disentangle the magnetic field properties, 
this will require other observables. 



3.1.2 Rotation measure 

Another independent approach to constrain magnetic field mod- 
els are Faraday rotation measurements. Theoretically, one expects 
the magnetic field in shocks to b e aligned with them due to shock 
compression (Enfili n etalJl 1998b and stretching and shearing mo - 
tions induced by oblique shocks l lSchekochihin & Cowlevll200d) . 
In combination with the small synchrotron emitting volume that 
is caused by the small synchrotron cooling time, this yields to 
polarized relic emission. Indeed, radio relics have been observed 
to be polarized up to the 40 per cent level iFeretti et al.l 120041 ; 
IClarke & E nfilin 2006). When polarized radio emission propagates 
through a magnetized medium, its plane of polarization rotates for a 



nonzero line-of-sight component of the magnetic field B- due to the 
birefringent property of the plasma - Faraday rotation. The Faraday 
rotation angle is given by 



A 2 RM + 0„ 



where 



RM(x ± ) = 



a I B : (x)n e (x)dz 
Jo 

rad B 



=a 812 - 



n. 



(9) 



(10) 



(11) 



m 2 11G 10 -3 cm Mpc ' 

where ao = e 3 /(27tm 2 c 4 ), x = (x x , I), and « c is the number density 
of electrons. In Eqn.QT] we have assumed constant values and a ho- 
mogeneous magnetic field along the line-of-sight to give an order 
of magnitude estimate for RM values. Assuming statistically homo- 
geneous and isotropic magnetic fields, the RM dispersion (RM 2 ) 
reads as follows, 



(RM 2 ) 



f 

Jo 



n e (x)B z (x)dz 



af^-L(n e (xfB z {xf) 
a 2 ^LC„ cBz (n e (x) 2 }(B z (x) 2 ), 



where 



C„ e B„ 



<gg) 
{B 2 ){nlY 



(12) 

(13) 
(14) 

(15) 



is defined as the correlation fact or and A B = 2/3 x A , is the 3D 
magnetic auto-correlation scale (Enfil in & Vogtl [2003) which can 
be estimated from the measured RM power spectrum. 

We studied the Faraday rotation of the largest relic in cluster 
g72a (Fig. [5}. The aim is to recover intrinsic statistical properties 
of the ICM magnetic field by studying RM statistics. We produce 
the RM map by projecting the line-of-sight component of our tan- 
gled magnetic field (Sect. l2~2t and the thermal electron density that 
has also been mapped from its Lagrangian distribution onto a 3D 
grid0 Firstly, we are concerned about the observability of polarized 
emission. The RM scales with By and n e according to Eqn.llOlsuch 
that we expect RM values to increase across the relic towards the 
projected cluster center. However, large RM values leads to con- 
fusion when tr ying to observe the polariza tion angle. Beam width 
depolarization ( Gard ner & Whiteoakl 19661) takes place if the polar- 
ization angle changes by a radian on scales shorter then the beam. 
To avoid this, one can go to shorter wavelengths and smaller beams 
at the expense of radio luminosity. Secondly, we are concerned with 
RM contamination from the galaxy. Even at moderately high galac- 
tic lat itudes the galactic RM contribution dSimard-Normandin et al.l 
1981) can be approximately the same order as the RM we calculate. 
However, the strength of our RM depends on the relic location with 
respect to the cluster and observer, as well as our magnetic field 
parameters. Therefore, a different location or parametrization will 
lead to stronger or weaker RM. Also, the galactic RM in principle 
can be modelled and removed from the RM map. 

High-quality rotation measure maps enable one to measure the 
RM power spectra. The large angular extent of giant radio relics 
provides a powerful tool of probing the maximum coherence scales 
of the magnetic field in clusters; in contrast RM maps from radio 



4 Note that we neglect gas above the very conservative threshold of number 
density n t hres = "sft/32 to be consistent with our proceeding in Sect. 12.4.21 
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Figure 5. Left panel: Faraday rotation measure map of the largest relic in g72a, if the cluster were placed at z~0.05 (mimicking A2256). Right panel: 
Polarization angle map at v = 1 .4 GHz, assuming a uniform rotation angle at the position of the relic. The cluster center is located in the direction of the upper 
left corner for both the images. Regions with high RM variance correspond to high spatial variation of the polarization angle. In combination with a finite 
beam size, this makes it challenging to observe a high degree of polarization in such a relic. The magnetic field parameters are org = 0.7 Bq = 2.5 \xG and the 
contours represent the surface brightness increasing in decades from 5 X 1CT 3 mJy arcmin~ 2 at v = 1.4 GHz. 



Table 4. Magnetic field realizations (MFR) used in RM maps. 



magnetic field 


correlation length 


input slope of 


realization name 


(A B ) [ft- 1 kpc] 


power spectra 


MFR1 


100 


-5/3 


MFR2 


200 


-5/3 


MFR3 


100 


-2 



lobes are typically much smaller. We calculate power spectra from 
our RM maps and magnetic field realizations for each model sepa- 
rately (Fig. [6] and Table |4j. For consistency reasons, we only con- 
sider the volume subtended by the radio relic when calculating the 
magnetic power spectrum. We define the RM power spectra (P[rm]) 
and the power spectrum of B- (P[b-]) as follows, 

(RM 2 ) = 2k f kP lRM] dk, (16) 
Jo 

(Er) = 3-4n f jfPwjdk, (17) 
Jo 

where the additional factor of three accounts for fluctuations in the 
total magnetic field while assuming our random Gaussian field. A 
partial Monte Carlo method was used to determine the lcr error 
bars on the power spectrum. Assuming a constant magnetic coher- 
ence scale, we construct the envelope function by computing the 
variance of RM (Eqn.[T3j- We multiply this envelope function with 
10 3 realizations of random Gaussian field and measure the power 
spectrum on each of these maps. The fractional errors are computed 
from the variance of these power spectra. While the first six power 
spectrum bins are on average 50 per cent correlated, the correla- 
tions drop to be below 20 per cent for the bins on smaller scales. 



By construction our parametrization of the magnetic field is 
correlated with the electron number density « e which might pos- 
sibly introduce biases in our RM maps. However, comparing RM 
power spectrum to the magnetic power spectrum, we find that the 
RM power spectrum overall resembles the original shape of the 
magnetic power spectrum and the injection scale corresponds to 
the scale of maximal power in the RM map. This is true for all our 
magnetic field realizations. We measure the slope of these power- 
laws at small scales to an accuracy of ±0.05, which is significant to 
differentiate between Kolmogorov (£~ 5 ' 3 ) and Burgers (k~ 2 ) turbu- 
lence spectra. We also find that the measurement of the RM power 
spectrum slope is independent of the magnetic field parameters a B 
and B . The full correlation matrix of the power spectrum bins 
is used when fitting the power-laws and in the error calculations. 
We find that our measured slope is flatter then expected and can 
be attributed to small scale fluctuations in n e , since P[ B ,] has the 
same slope as the power spectrum of B- integrated along the line of 
sight. The small inhomogeneity of n e in our simulations does not 
severely affect the intrinsic spectral shape. Thus, it is possible in 
principle to recover the intri nsic 3D magnetic power spectrum by 
solving the inverse problem dVogt & EnBl"inll2003l : lEnfilin & Vogtl 
l2003l : IVogt & Enl31inll2005h . 

Using Eqns. [14] Q6] and [T7] we estimated the rms magnetic 
field strength from P[rm\ and Py B -\ respectively and recover our ini- 
tial rms magnetic field strength. We find that our correlation factor 
( s/C^b- — 5.3)0 is 20 per cent larger than the correlation factor 



5 We caution the reader that the particular value of C„ c b- reflects the 
parametrization of the magnetic field we adopt in our model, and may be 
realized differently in Nature. We also note th at the value of the cor re- 
lation factor in the non-radiative simulation by Pfrommer et al. 12008) is 
^jC ni! B- — 6- Further work is required to address this question in the context 
of MHD cluster simulations. 
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Figure 6. Top: power spectrum of the RM map P[ R M](k) and power spec- 
trum of the line-of-sight component of the magnetic field P[B.](k) for 
MFR1, where the error bars represent the Icr confidence regions. P{RM](k) 
attains excess power at large angular scales from fluctuations in n e . Our 
RM power spectrum matches the shape and the peak scale of the input 
power spectrum within the error bars. Bottom: power spectra of RM maps 
for different magnetic field realizations (cf. Table[4}. All RM power spec- 
tra recover the shape and characteristic scale of their magnetic input power 
spectra. 



obtained by fitting a smooth y3-model to the spherically averaged 



profile of n c and scaling B • 



with the same a B = 0.7 that we 



used to construct our RM maps (similar to the procedure applied by 
lEnfflin & Vogtll2003l and lMurgia et al.ll2004h . This result suggests 
that the fairly homogeneous density distribution in our simulations 
(after removing the galaxy contamination described in Sect. 12.4.21 
does not severely bias the average magnetic field strengths esti- 
mated by RM studies if one takes into account the overall shape of 
the profiles of n c and B. 

For convenience, we derive a formula for the rms magnetic 
field strength ( -\/{B 2 )) as a function of the peak of P[ RM ] and the 
rms fluctuations of « c . Multiplying P[rm\ with a Heaviside function 
&(k) and ensuring that the intrinsic spectrum is sufficiently steeper 
than lc 2 , Eqn.ll6lbecomes 



{RM 1 ) 



kP[RM] (k) &(k- fcpeak) dk, 



In 

nP[RM\ (Apeak) 



(18) 



The value of P[£M] (Apeak) A 2 cak can be read off directly from Fig. [6] 
and combining Eqns. [T5] and [14] yields an approximate value for 



2nPy RM] (fcpeak) A pciik 

alA B LC„ tB , (nj) 



O.lfiG 



' P[RM](k p ^) k 



2 

peak 



50radm 2 



10" 4 h 2 cm- 3 j 
i i 



(20) 



,100/r'kpc/ \4/7-'Mpc/ \ 27 

where we inserted numerical values from our simulation in the last 
step. 

3.2 Existence and properties of the WHIM 

In this section, we investigate the potential of radio relic observa- 
tions to infer the hydrodynamic properties such as density and tem- 
perature of the WHIM. 



3.2.1 Properties of virializing shocks 

Diffusive shock acceleration determines the shape of the CRe spec- 
trum that we model as a power-law momentum spectrum (neglect- 
ing non-linear effects). Synchrotron losses cause a steepening of 
this power-law by one power of momentum. Spatially inhomo- 
geneous virializing shocks with a distribution of shock strengths 
cause a spatial variation of the spectral index of the cooled CR elec- 
tron spectrum. This is reflected in an inhomogeneous distribution of 
synchrotron spectral index that may help to reconstruct the merg- 
ing geometry by providing a snapshot of the structure formation 
process in a galaxy cluster. 

The spectral indices of the radio surface brightness a v yx> and 
that of the intrinsic 3D emissivity a v ,3D are defined by, 



Q-V.2D = - 



OtySD 



log(T^) 
" log(f) ' 



(21) 



(22) 



(19) 



It is unclear ab initio whether the projected spectral index repre- 
sents the actual deprojected quantity (0^30) due to possible su- 
perposition of different radio emitting structures along the line-of- 
sight. We study how these two quantities relate to each other, and 
present 2D spectral index maps of both the largest radio relic in 
g72a and the total emission from the same area (cf. Fig. |7J- We 
note that more than 99 per cent of the total radio emission can be 
attributed to emission within the radio relic. As a result, the a Vt 2D 
maps are not contaminated by the diffuse spurious emission. Note 
the edge effects in Fig.|7]that show up in the projection of the sin- 
gle relic where the emission falls off. These effects are due to the 
sharp emissivity cutoff of our relic finder and incomplete sampling 
of SPH relic particles at different frequencies. In regions with high 
synchrotron brightness, one can ignore these edge effects, and the 
resulting distribution of or v2D is fairly uniform ((o^d) - 1.15, with 
<t »v2d ~ 0.04) implying that this relic traces a single structure for- 
mation shock wave. 

We further study the distribution of or vZD and Q' V 3 D in our 
largest individual relic to uncover a connection between them. 
Probability distribution functions (PDF) are constructed for both 
*v,2D and «v,3D for varying parameters of the magnetic field (cf. 
Fig-H}- To avoid contamination from edge effects seen in Fig.|7] the 
ffy,2D PDF was made for the brightest 20 per cent of the pixels and 
the a v3D PDF was weighted by particle emissivity. These distribu- 
tions do not change with our choices for magnetic field parameters 
implying that the spectral indices are practically independent of the 
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Figure 7. Spectral index map, a v 2D, between 150 MHz and 1.4 GHz for the largest relic in g72a, with only the SPH particles selected by the relic finder 
projected (left panel) and for the total emission in the same region (right panel). The contours show orders of magnitude in surface brightness in mJy arcmin~ 2 , 
with the highest contour representing 5 mJy arcmin~ 2 at y = 1 .4 GHz. Notice the edge effects that show up in the projection of the single relic where the 
emission falls off. These effects are due to the sharp emissivity cutoff of our relic finder and incomplete sampling of SPH relic particles at different frequencies. 
More importantly, in regions with high synchrotron brightness, the spectral index is almost uniform across the central relic implying that this relic traces a 
single formation shock wave. 



magnetic field and depends mainly on properties of the shock. An- 
other striking result is that the median values for cr^D and or^o are 
statistically consistent within 1-cr. Assuming that the line-of-sight 
integral is dominated by one bright relic and choosing a pixel scale 
that is smaller than the length scale on which the post-shock den- 
sity varies, we can easily show that the 2D and intrinsic 3D spec- 
tral index are identical. If there are more radio emitting regions 
contributing to the observed surface brightness, we expect a con- 
cave radio spectrum. Synchrotron cooling as well as re-acceleration 
lead to spectral steepeni ng in particular at high radio frequencies 
dSchlickeiser et alj|l987h . Future work is required to address the 
associated biases of the relation between a v iu and a viD . 

Figure [9] shows the observable parameters space of relic lu- 
minosity and the median 3D spectral index. This parameter space 
compares the shocks strength, which is related to the 3D spectral 
index (cf. Eqn. IClt to the energy dissipated at the shock, which is 
related to the relic luminosity. There is a trend that strong shocks 
are associated with the more luminous relics. The implications of 
this trend is that the brightest radio relics should show predom- 
inantly flatter spectral ind ices, which is the current observational 
status of giant radio relics dFerrari et alj|2008h . 



3.2.2 Predicting pre-shock properties 

A majority of the hot gas (> 10 7 K) found at the centers of galaxy 
clusters is believed to originate from the WHIM that is shock heated 
through large-scale virializing structure formation shocks. These 
structure formation shocks are traced by synchrotron emission in 
form of radio relics from recently accelerated electrons (cf. Sect. 
|2.3l l. We have shown that under particular conditions, the observed 
median 2D spectral index corresponds to the weighted median 3D 
spectral index a F ,3D (cf. Sect. l3.2.Tt . The 3D spectral index can be 
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Figure 9. Two dimensional observable parameter space for radio relics in 
galaxy cluster g72a. Each symbol represents a relic within g72a and they 
are characterized by total luminosity and median of the 3D spectral index. 
Shown are observable relics (red crosses), theoretical relics (blue circles) 
and the large relic is emphasized by the bold cross and the filled circle. 
There is a trend to high spectral indices and a larger scatter for less lumi- 
nous relics. As shown in Fig. [4] the loss of the less luminous relics from the 
observable to theoretical parameters is a result of the brightest relic swal- 
lowing up these smaller relics due to a decrease of the emission threshold 
parameter. 



related to the Mach number of the shock (cf. Eqn. IClt under the 
assumption that we have an ideal fluid with a given adiabatic in- 
dex. One can obtain information on the post-shock values for den- 
sity, pressure and temperature of the ICM through deprojections 
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Figure 8. Top: spectral index distributions for our largest relic using our standard particular magnetic field parametrization. 2D spectral index distributions, 
f(@v2D), of the full map and the brightest 20 percent pixels (top left) are contrasted to the radio luminosity weighted distribution of the 3D spectral index, 
dJ va /da v 3E) (top right). Bottom left: the median of f(a v 2D) for the brightest 20 percent pixels (S v > 12 mjy arcmin~ 2 at v = 150 MHz and S v > 1 mjy 
arcmin -2 at v = 1.4 GHz) as a function of the magnetic decline ag with the error bars representing the lcr percentiles. Bottom right: the median of dj va /da!y,3D 
as a function of the magnetic decline ag with the error bars representing the lcr percentiles. This shows that for the giant radio relic of g72a, the median of the 
2D and the 3D spectral indices agree statistically and are almost independent of the magnetic decline ag- 



of de ep X-ray or Sunyaev-Zel'dovich observations dZaroubi et al.l 
1998). With the knowledge of the Mach numbers combined with 
post-shock values we calculate the pre-shock conditions of the ICM 
(the WHIM) using the Rankine-Hugoniot jump conditions (cf. Ap- 
pendixjCt. 

We take an optimistic approach and assume that the deprojec- 
tions of the thermal observables can be done ideally such that we 
use our radial profiles calculated from the solid angle subtended by 
the largest relic for simplicity (cf. Fig. 1 10t. We define the shock re- 
gion by locating the radial bins that contain the majority of shocked 
relic particles (> 85%). A small fraction of the relic particles leak 
into radial bins adjacent to the shock causing slight enhanced values 
of the radial profile. As mentioned above we insert the calculated 
average Mach number and the post-shock value into the Rankine- 
Hugoniot jump conditions for density, pressure and temperature to 
estimate the upper limits of these WHIM properties. Our predicted 
upper limits for pressure and temperature of the WHIM are con- 



sistent with the simulated pre-shock properties within one standard 
deviation. We note that the standard deviation of these hydrody- 
namic properties reflect actual physical variations due to an oblique 
shock that is not perfect tangential. Additionally, the particular relic 
chosen is located at ~ R 2 m/2, which is within the cluster volume. 
There are other observat ionally know relics that reside at the virial 
radius and beyond (e.g. iBagchi et alj|200rj) . These relic are better 
suited to probe the WHIM in combination with future X-ray and 
multi-frequency SZ data. Thus, this example is to be taken as a 
demonstration of our concept. 

In the following we want to address possible biases with our 
method and show that the discrepancy between the predicted val- 
ues and the average radial value of the WHIM can be explained 
by differences between the calculated Mach numbers and the me- 
dian of the weighted Mach numbers in the radio relics (cf. Fig. 
lilt . We find that the weighted Mach numbers have systematically 
lower values compared to the theoretical expectation due to the 
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Figure 10. Radial profiles of galaxy cluster g72a restricted to the solid angle 
subtended by the largest relic for the density, pressure and temperature. The 
shocked region is seen in the profile at roughly R200/2 and is marked by 
the black diamonds with one sigma error bars, except for the density profile 
where the error bars are too small to show. The predicted pre-shock values 
(crosses) are ~ 7%, ~ 67% and ~ 67% different from the average profile 
values, but fall within the standard deviation for pressure and temperature. 
Most of this variation is caused by the shock being oblique and not perfectly 
tangential. 



skewed distribution of the emissivity weighted a v3D . According to 
the Rankine-Hugoniot jump conditions, systematically higher val- 
ues of the shock strength should over-estimate the jumps and hence 
under-predict all the pre-shock quantities, which appears to be the 
case (Fig.llOt, 

3.3 Dependence on dynamical state and cluster mass 

We study the distribution of radio relics for the entire galaxy clus- 
ter sample, which shows a variety of both dynamical states (ranging 
from merging to cool core clusters) and masses (a range of almost 
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Figure 11. Top: median of the Mach number distribution for each relic as 
a function of the median of the distribution function of ovjd- Bottom: the 
observationally relevant quantities are the weighted distribution functions. 
Shown are median of dejiss/dAI as a function of d./ vn /dav,3D, where each 
point represents a radio relic in cluster g72a. Shown are observable relics 
(red crosses), theoretical relics (blue circles) and the large relic is empha- 
sized by the bold cross and the filled circle. The dashed line is the theo- 
retical Mach number computed from directly from the av,3D (cf. Appendix 
Eqn. ICU . The theoretical Mach number over-predicts the median of the 
weighted Mach number, due to the weighting of a y ^o by the skewed distri- 
bution function of the radio luminosity. 



two orders of magnitude). In Fig.[T2] we investigate how cluster 
mass and dynamical state depend on the relic luminosity function. 
In the case of our theoretical parameter space, more massive galaxy 
clusters clearly have more radio relics than the lower mass clusters 
with a power-law scaling of M^£. (Fig. 1 13t. For the current obser- 
vational capabilities, we predict that only the most massive clusters 
should have a significant sample of radio relics. 

Ideally, one would like to directly compare clusters with the 
purpose of using relic number statistics as a mass proxy. However, 
the luminosity functions have another trend, which causes the scat- 
ter in the relationship between clusters mass and total number of 
relics per cluster. This trend relates the cluster's dynamical state to 
the luminosity of their brighter relics. The clusters g72a and g51 
have the same virial mass 1.1 x 10 15 Mq, but g51 is a relaxed cool 
core clusters in contrast to the active state of g72a. One can see the 
two most luminous radio relics of g72a are an order of magnitude 
brighter than any of g51's relics. Furthermore, the total amount of 
relics of g72a is greater than that of g51. Merging clusters inher- 
ently have more shocks which yields to more high-energy CRes 
and magnetic field amplification resulting in more radio relics. This 
trend is even more severe for galaxy clusters of smaller mass (Fig 
1 1 2b - Our results show a larger probability of observing a relic in a 
more massive cluster that is dynamically active. This dependence 
on mass and dynamical state offers a possible explanation for why 
all current observed radio relics are in massive merging clusters. 
They are expected to be the brightest of a dimmer population of 
radio relics. 
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Figure 12. Luminosity functions for our sample of 10 clusters with the high mass clusters (M v j, > 2 X 10 hT Mq) in the upper panels and the low mass 
clusters in the lower panels for our standard magnetic model at 150 MHz (see Table[2}. The left panels contain relics found using the observable parameters 
and the right panels contain the theoretical parameters. The luminosity functions from the observable parameters show that more massive clusters have relics 
while low mass clusters have no relics (with gib and g72b being the exceptions in both cases, respectively). The cumulative luminosity functions from the 
theoretical parameters show the trend for higher mass clusters to contain more relics and dynamical cluster stage modulates this effect notably, especially at 
low cluster masses. For instance compare the luminosity functions of the equal mass clusters g51 and g72a, the latter of which is a post-merging cluster. 



4 DISCUSSION 

4.1 Comparison with previous theoretical work 

Previously, there has been analytical work feshet et alj2004ah and 
pioneering cosmological simulations lMinia tietalj200lL 2000) on 
studying cluster synchrotron emission from shock-accelerated elec- 
trons. The latter authors simulated the non-thermal cluster emis- 
sion by numerically modelling discretised cosmic ray (CR) en- 
ergy spectra on top of Eulerian grid-based cosmological simu- 
lations. Their approach neglected the hydrodynamic pressure of 
the CR proton component, was quite limited in its adaptive reso- 
lution capability, and neglected dissipative gas physics including 
radiative cooling, star formation, and supernova feedback. To al- 
low studies of the dynamical effects of CR protons in radiatively 
cooling galactic and cluster environments, a CR proton formalism 
was developed that is based on smoothed p article hydrodynamical 
representation of the equations of motion jPfrommer et al J [20061 ; 



lEnfflin et al.ll2007l : ljubelgas et aHl2008h . The emphasis is given to 
the dynamical impact of CR protons on hydrodynamics, while si- 
multaneously allowing for the important C R proton injection and 
loss processes in a cosmological setting. IPfrommer etai] (2008) 
modelled the CR electron components due to shock acceleration 
as well as those being produced in hadronic CR proton interac- 
tion w ith ambient gas protons. Us ing t his formalism an d mod- 
elling, IPfrommer et ail < l200l 120081) and lPfromm er (2008) coher- 
ently studied the non-thermal cluster emission, the cosmic ray pres- 
sure component, and its implications for thermal cluster observ- 
ables such as the X-ray emission and the Sunyaev-Zel'dovich ef- 
fect. The focus of iMiniati et al. I d200lh was on primary CRe syn- 
chrotron emission from galaxy clusters as a whole, whereas we fo- 
cused in this work on the emission from individual relics in detail 
to study the synchrotron observables and how they are sensitive to 
the large scale magnetic properties. To summarize, we have im- 
proved on past work in the observable predictions, simulations and 
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Figure 13. Relation between the number of radio relics within a galaxy 
cluster and the cluster mass, where each point represents one of our sim- 
ulated clusters and the line represent the best fit power-law that scales as 
A/?™. There is a trend for higher mass clusters to contain more relics. 



observational understanding of non-thermal emission from primary 
accelerated CRe. 

Additionally, we point out that IC and y-ray emission are al- 
ternate ways to study struct ure formation shocks. T his was pro- 
posed in analytical work by lLoeb & Waxmar i feOOtt) and in sim- 
ulations dMiniatil 120031: iKeshet et al.1 120031 : iMiniati et al.1 120071 ; 
IPfrommer et al.l2008l ; |Pfrommej2008h . We are optimistic that high 
energy y-ray experiments, such as the Fermi y-ray space telescope 
(formely GLAST) and future imaging air Cerenkov telescopes, will 
aid in further developing the picture of non-thermal emission at 
structure formation shocks. 



4.2 Assumptions 

In our attempt to model diffuse radio relic emission from galaxy 
clusters, we have made several simplifying assumptions. (1) We as- 
sume the modified thermal leakage model fully describes the pro- 
cess of diffusive shock acceleration and did not vary the param- 
eters associated with it. The described observations allow one to 
test the self-consistency of this hypothesis, and might finally allow 
improvement of our knowledge about diffusive shock acceleration 
in high-/? plasmas. (2) We neglect at this point the modifications 
of this simple model due to non-linear shock acceleration, as well 
as cosmic ray modified shocks, and postpone their study until fu- 
ture work. (3) We also neglect re-acceleration of mildly relativis- 
tic electrons that have been injected in the past either by forma- 
tion shocks or other sources such as AGN. (4) We use a simple 
parametrization for the magnetic field. There are indications that 
the m ain characteristics of this model are realiz ed in clusters on av- 
erage iVogt & EnBlinl20 05 ; Murgi a et al.l2004h . Future work has to 
be dedicated to study the distribution of magnetic fields that follow 
the magneto-hydrodynamics in radiative simulations. (5) We are 
solving for a steady-state spectrum of the electron population and 
are not sensitive to spectral aging processes across the relic as they 
may have been found recently bv lGiacintucci et ail J2008h - (6) In 
our model, we assume the thermal reservoir to be the source of elec- 
trons. (7) In our analysis, we only consider the rotation measure sig- 
nal from the line-of-sight integration of the density weighted paral- 
lel magnetic field. We explicitly neglect possible contributions from 



magnetic field a mplifications due to post-shock turbulence local to 
the shock wave dVladimirov et alj|2004 pRvu et al.ll2008l) . These 
questions are beyond the scope of this work and will be studied 
elsewhere. (8) The mass contained within the relics is not a physi- 
cally relevant quantity and suffers from the finite resolution of the 
SPH technique at the dilute shocks in the virial regions of clusters 
and beyond. The radio luminosity, however, is a robust prediction 
within a given acceleration model since it reflects conserved quan- 
tities such as energy and mass across the shock. 



5 CONCLUSIONS 

The intermittency and inhomogeneous nature of structure forma- 
tion shocks are characterized by a highly non-Gaussian distribution 
function. This requires numerical simulations to study the implied 
non-universality of the induced radio relic (or gischt) emission. It is 
hard to conceive of an adequate analytical approach to this problem. 
Observing the polarized emission of a sample of relics at differ- 
ent frequencies enables us to gain insight into the non-equilibrium 
processes at work - in particular into the interplay of large scale 
magnetic fields and structure formation shocks. The relevant ob- 
servables of the relics include morphology, spectral shape, relic lu- 
minosity function and Faraday rotation measure. The theoretical 
implications of radio relic observables are as follows: 

We model the shock acceleration of electrons at formation 
shocks and find that the morphology of radio relics unambiguously 
characterises the underlying structure of dissipating shock waves 
(Fig. 2). The resulting simulated relics are very similar to the ob- 
served relics and thus support our hypothesis. Their positions iden- 
tify regions that are not in equilibrium and where the electron and 
ion temperatures are expected to show strong deviations due to 
the comparatively long Coulomb mean free path that governs their 
equilibration process. 

The relic luminosity function is sensitive to the combination 
of normalization and scaling properties of the magnetic field with 
thermal energy density as well as the electron acceleration effi- 
ciency, the mass and dynamical state of a cluster. Thus, it can pro- 
vide hints about the processes that generate these large scale fields 
and can help to disentangle the dominant transport processes which 
include effects from magnetic flux freezing and growth by turbulent 
dynamos. 

The rotation measure (RM) map is sensitive to the line-of- 
sight integrated magnetic field. From the mean and variance of RM 
maps, we can infer the location of the formation shock with respect 
to the cluster center as the variance increases as a function of inte- 
gration length (Eqn.ll4t. This helps in constraining the geometry of 
the merger. Deprojecting the 2D RM power spectrum enables one 
to measure the 3D magnetic power spectrum, under the assumption 
that the behaviour of the electron density along the line of sight can 
be obtained from X-ray measurements. The peak of the 3D power 
spectrum yields the total magnetic energy and the magnetic coher- 
ence length A B . Performing this procedure for different relics or for 
different regions of one large relic allows us to estimate the vari- 
ance of A B across the cluster and might possibly tell us about the 
nature of MHD turbulence. We found that the correlation between 
« c and B biases the rms magnetic field strength derived from RM 
maps high if this is not taken into account; we note that most works 
have done so. This correlation should be a natural consequence of 
MHD effects such as flux freezing. If systematic errors associated 
with RM studies are smaller than statistical ones, we find that mea- 
surements of the small scale slope of the RM power spectrum are 
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accurate enough to differentiate between Kolmogorov (k~ 513 ) and 
Burgers (Ar 2 ) turbulence spectra. The interpretation of these slopes 
is however not straight forward and needs to account for an addi- 
tional flattening due to small scale fluctuations in n c . 

The spectral index of a power law spectrum of a radio relic is 
a measure of the shock strength of that relic. The median spectral 
index of a distribution of relics probes a distribution of virializing 
shocks and can give an indication of CR proton injection. This is of 
particular relevance for questions concerning the pressure contribu- 
tion of non-thermal components and enables comparison with pre- 
dictions of hydrodynamical simulations. The shape of the spectrum 
is sensitive to the acceleration mechanism of the relativistic elec- 
trons and to their cooling processes. The variation of the spectrum 
over the relic allows one to infer in situ magnetic field strengths by 
comparing the synchrotron and IC cooling times to the advection 
time downstream provided that the magnetic energy density is not 
much smaller than the CMB energy density. This might constrain 
models for the magnetic amplification at shocks in high beta plas- 
mas. A radio relic's luminosity is roughly correlated with the shock 
strength. Thus, it is favourable to look for radio relics in large, dy- 
namically disturbed clusters, or use relic detections as a proxy for 
dynamical activity of clusters (Schuecker et al. 2 00 lh . 

We demonstrated that the combination of the relic spectral in- 
dex with deprojected X-ray and SZ profiles allows one to indirectly 
infer upper limits on the density and temperature of the warm hot 
intergalactic medium. Simulations show that the WHIM is not uni- 
form, rather it is characterized by highly inhomogeneous structure 
that shows intermittent accretion events that are channeled mostly 
through filaments. 

We predict that there will be a large sample of polarized radio 
relics from a considerable number clusters in the near future. This 
sample should allow one to constrain macroscopic model param- 
eters, which are expected to be higly non-Gaussian, using a joint 
analysis method on the radio observables. For example, the com- 
bination of the relic luminosity function, RM power spectra, X- 
ray and SZ measurements should constrain the models of diffusive 
shock acceleration and large scale magnetic fields. Future work will 
address the details of this procedure. 
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APPENDIX A: INTERPOLATING AND PROJECTING 
SPH QUANTITIES 

In the course of this work we are required to interpolate our La- 
grangian energy density distribution as given by SPH on a 3D grid. 
We remind the reader that the SPH smoothin g kernel of an SPH par - 
ticle a, W(|r-r„|, h a ), is given by Eqn. A.l o flSpringel et al.l ( f200lh . 
It is normalized in the continuum such that f W(r,h)d 3 r = 1. A 
scalar field x(rf^ is interpolated onto a 3D grid cell at ryj by the 
product of itself with the specific volume M a jp a of the gas parti- 
cles over a comoving cube, 



(Al) 



where L 3 . is the comoving volume of the grid cell and we define 
the normalized 3D smoothing kernel of SPH particle a at the grid 
position r iik by 



<W(\r m -r a \,h a ): 



W(\r, jk -r a \,h a ) 



(A2) 



Yjijk W(\r ijk -r a \,h a ) 

We note that the normalized interpolation conserves the interpo- 
lated quantity strictly without any further requirement on the grid 
size. 

Similarly, we employ the method of normalized projection of 
a three dimensional SPH scalar fields x(r) to perform projection 
integrals yielding the quantity X(r ± ). In analogy to Eqn. IA1I we 
obtain 



L 2 . ■^- J p„ 

nix n ' u 



■r a \,h a ), 



(A3) 



where Z? ix is the comoving area of the pixel and the normalized 2D 
projected smoothing kernel of SPH particle a at the grid position 
r ±J j derives from the projected SPH kernel Fflr^yl, h„) and is given 
by 

. . . Y(\r ± ij — r a \, h a ) 

y(\r±, u - r„\, K) = 1 — . (A4) 

Zij k y(\r±jj -r a \,h a ) 



APPENDIX B: THEORETICAL EMISSION THRESHOLD 

The emission threshold for the observable parameters and theo- 
retical parameters differ by 12 orders of magnitude. This dynamic 
range is beyond the ability of any future telescope on the horizon. 
Varying the emission threshold of our theoretical parameters by 
six orders of magnitudes only very weakly affects our results. In 
particular, we show in Fig. IB II that such a dramatic variation has 
only little influence on the high-end of the radio relic luminosity 
function. 



6 We note that in general, x has to be a thermodynamic extensive volume 
density such that the product xM/p is extensive. 
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Figure Bl. Cumulative luminosity functions for different values of the 
emission cutoff. This shows the robustness of our predictions for future in- 
strument capabilities. 



APPENDIX C: RANKINE-HUGONIOT CONDITIONS 

The three dimensio nal spectral index c an be transformed into a 
Mach number (M) tenfllin et al.ll2007l) . if one assumes an ideal 
fluid that is characterized a a single adiabatic index y, 



4(l + a v , 3D ) 
Al l+4a v , 3D -3 r 

Under these conditions, the well-known Rankine-Hugoniot jump 
conditions allow to relate the hydrodynamic post-shock quantities 
(denoted with a subscript 2) to the pre-shock quantities (denoted 
with a subscript 1), 

pi (y-l)M 2 + 2 

p 2 (y + l)M 2 ' 1 ; 

r, (y + l) 2 M 2 + 2 



T 2 " [2yM 2 -(y-l)][(y-l)M*+2Y 
Pi 7+1 



(C3) 
(C4) 



P 2 2yM 2 -{y-\) 

Phenomenologically, we show in Fig. QT] that Eqn. IC1I under- 
predicts the average {M) if one were to infer M from spectral in- 
dex maps. This translates into an upper limit for the predicted pre- 
shock density. In the case of temperature and pressure, the under- 
prediction of the average (M) leads to an over-estimation of the 
pre-shock values which translates into lower limits for temperature 
and pressure (Eqns. lC3l and lC4b . 
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